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Abstract 

We investigate the detectability of neutralino Dark Matter via direct and indirect 
searches as well as collider signatures of an 50(10) model with two intermediate scales. 
We compare the direct Dark Matter detection cross section and the muon flux due to 
neutralino annihilation in the Sun that we obtain in this model with mSUGRA predic- 
tions and with the sensitivity of current and future experiments. In both cases, we find 
that the detectability improves as the model deviates more from mSUGRA. In order 
to study collider signatures, we choose two benchmark points that represent the main 
phenomenological features of the model: a lower value of and reduced third genera- 
tion sfermion masses due to extra Yukawa coupling contributions in the Renormalization 
Group Equations, and increased first and second generation slepton masses due to new 
gaugino loop contributions. We show that measurements at the LHC can distinguish this 
model from mSUGRA in both cases, by counting events containing leptonically decaying 

bosons, heavy neutral Higgs bosons, or like-sign lepton pairs. 
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1 Introduction 



Grand Unified Theories (GUTs) based on the gauge group 5*0(10) have been considered good 
candidates for the unification of electroweak and strong interactions [T]. All matter fields of 
one generation are incorporated in a single irreducible representation, the spinor 16. Moreover, 
the "seesaw" mechanism [2] , which can explain small neutrino masses as indicated by neutrino 
oscillations, is naturally embedded. 

In [3], two of us in particular chose a model by Aulakh et al. a supersymmetric 50(10) 
model with two intermediate scales: 50(10) is first broken to SU{4:)c x SU{2)l x SU{2)ji by 
a 54 dimensional Higgs at the GUT scale M^; then to 5?7(3)c x U{1)b-l x SU{2)l x SU{2)r 
by 45 at scale Mc] finally to the Standard Model gauge group by 126 + 126 at scale Mr. 
Imposing the unification condition for the gauge couplings fixes the intermediate scales Mc 
and Mr for given Mx] i.e. Mx is a free parameter. However, its lower bound is set by 
the lower bound on the lifetime of the proton [5]. We took Mx = 3 ■ lO^^GeV as default 
value. In addition, a second pair of Higgs doublets was introduced, in order to modify the 
minimal 5*0(10) predictions for the masses of quarks and leptons, which are not consistent 
with experiments. In order to compare the low energy phenomenology of the model with that 
of mSUGRA [6], we assumed universal soft breaking parameters (mo,Mi/2,^o) as boundary 
condition at the GUT scale. 

The 126-dimensional Higgs whose vacuum expectation value breaks SU{2)r x U{1)b-l 
to U{1)y also gives Majorana masses to the right-handed neutrinos. The resulting masses for 
the light neutrinos are schematically written as 

^ = ^ = iMH^' a) 

" Mn Ym{12Q) ■ ^ ' 

The Yukawa coupling Yat, and to a lesser extent Y^, gives new contributions to the Renormal- 
ization Group Equations (RGEs) of the MSSM Yukawa couplings and soft breaking parame- 
ters. Therefore the weak-scale masses, and thus the radiative electroweak symmetry breaking 
and the relic density of Dark Matter, depend on Yjv, and hence on the light neutrino masses 
for fixed Mr. This remains qualitatively true for other 5*0(10) GUTs with a "type-I" seesaw 
mechanism at an intermediate scale. Note that unifies with the up-type quark Yukawa 
couplings, and is hence fixed. For given Mr, and hence given (126), the absolute scale of the 
light neutrino masses is thus determined by Yat, with larger Y^ yielding %/iter neutrinos. For 
our minimal choice Mx = 3 • 10^^ GeV, the requirement that Y^ remains perturbative at least 
up to scale Mx therefore leads to the lower bound my > 0.15 eV on the mass of the heaviest 
light neutrino. 

For our current study the most important modification of the weak-scale spectrum is the 
reduction of the higgsino mass parameter which comes about as follows. The weak-scale 
stop masses are reduced compared to the mSUGRA prediction, due to the Yukawa coupling 
given by [3] 

M^Yuk,422 3 \Yx {F'^^rF' + FElF) . (2) 

Here PVYuk,422 is the superpotential valid between the scales Mc and Mx, and S/, are 
in the (10,1,3) and (10,3,1) representation, respectively, of the gauge group SU{4:)c x 
SU{2)i X SU{2)r, and F and F'^ denote quark and lepton superfields in the (4,2,1) and 
(4,1,2) representation. This reduces the term oc Y^^ in the RGE for m'j^^, leading to an 



increase of m|^^ at the weak scale. As a result, electroweak symmetry breaking requires 
sma//er values than in mSUGRA. 

Another distinctive feature of the model is the rapid increase of the gauge couplings at high 
energies. This is due to the introduction of large additional Higgs representations, needed in 
order to break the gauge symmetry. As a result, relations between weak-scale and the GUT- 
scale soft breaking parameters are modified [3]. In particular, for a given universal gaugino 
mass Ml/2 at the GUT scale, the SO{10) model predicts much smaller gaugino masses at the 
weak scale than mSUGRA does. 

In ref.[3] it was shown that thermal neutralino Dark Matter remains viable, although the 
allowed region of parameter space is even more highly constrained than in mSUGRA. In the 
next Section, we will focus on the detectability of this Dark Matter candidate by direct and 
indirect searches. In Sec. 3, we compare signatures at the LHC between mSUGRA and this 
50(10) model for two benchmark points. Finally, we conclude in Sec. 4. 

2 Direct and indirect detection of Dark Matter 

All results presented in this Section are obtained using a modified version [3j of SOFTSUSY 2 . 
[7] to evaluate the mass spectra at the weak scale. These are then fed into micrOMEGAs 2.2 
[H] to calculate the LSP relic density. If this is found acceptable, we feed the same low-energy 
spectrum into DarkSUSY 5.0.2 [9j for the calculation of various Dark Matter detection rates 

m- 

Specifically, we compute: elastic LSP-proton scattering cross sections due to spin-inde- 
pendent as well as spin-dependent interactions; the muon flux resulting from LSP annihilation 
in the Sun; and the antiproton flux from LSP annihilation in the halo of our galaxy. In all 
cases we compare with the sensitivities of the best current and/or near-future experiments. 
In most cases, we only consider parameter sets leading to a relic density within two standard 
deviations of the value found by combining WMAP data with other cosmological observations 

VLcDuh'^ = 0.1131 ± 0.0068 {2a range). (3) 

We begin our discussion by analyzing the impact of reducing Mx, i.e. "switching on" 
the intermediate scales, on the elastic LSP-proton scattering cross section. Figure [1] shows 
the spin-independent (SI; left frame) and spin-dependent (SD; right frame) contributions to 
this cross sections as we vary Mx while most soft breaking parameters as well as the light 
neutrino masses are kept fixed. The scalar mass mo is varied along with Mx such that the 
relic density lies in the range of Eq.([3]). Note that mg ^ in this plot, i.e. we are in the 
region of significant higgsino-neutralino mixing, which is most favorable for direct neutralino 
Dark Matter searches. As a result, the spin-independent cross section is always well above 
the projected sensitivity of the XENONIOO experiment [12], while the spin-dependent cross 
section lies above the projected sensitivity of the DMTPC experiment |12j . 

While these gross features remain unchanged, we see that, as Mx decreases, i.e. as the 
deviation from mSUGRA becomes larger, the cross section is enhanced, so that for the low- 
est Mx, the spin-independent cross section slightly exceeds the limit set by the CDMS II 
experiment [13]. This is mostly due to the reduction of the LSP mass for fixed M1/2 in the 
model with intermediate scales, which we mentioned near the end of Sec. 1. In particular, for 
Mx < 8 ■ 10^^ GeV, m^o < rrit, so that XiXi ~^ tt annihilation is forbidden. The loss of this 
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Figure 1: The Mx dependence of the SI (left frame) and SD (right frame) neutrahno-proton 
scattering cross section, mo is varied such that the rehc density satisfies the constraint ([3]). The 
neutrahno mass = 124 GeV, for Mx = 10^^-^ GeV, increases to 235 GeV for Mx = lO^^-^ 
GeV. The experimental constraints are taken for m^o = 120 GeV. 



important annihilation channel has to be compensated by increasing bino-higgsino mixing, 
i.e. by decreasing fi, which in turn is accomplished by increasing ttiq. This leads to increased 
couplings of the lightest neutralino to neutral Higgs bosons as well as to the boson. Note 
that in scenarios with gaugino mass unification, first generation squarks are always much 
heavier than the lightest neutralino, suppressing their contributions to LSP-nucleon scatter- 
ing. Choosing mg 3> ^1/25 done here, further strengthens this hierarchy, so that Higgs and 
Z° exchange contributions largely determine the SI and SD cross sections, respectively. 

The curves in Figs. [1] show a noticeable negative slope even away from this threshold. In 
case of the SI cross section, this is due to the reduction of the mass of the heavier neutral CP- 
even Higgs boson with decreasing Mx, which goes along with the reduction of the weak-scale 
gaugino masses (although for fixed mo the ratio m^/m^o slightly increases with decreasing 
Mx i3j). Note also that decreasing m^o requires a simultaneous, if slower, decrease of /i, 
since otherwise the higgsino-component of Xi would become too small, yielding too small an 
annihilation cross section. This decrease of both weak-scale gaugino masses and of fi with 
decreasing Mx implies that the higgsino components of the LSP become more different in 
magnitude; note that they become identical in size for |/i| ^ Mz- This in turn enhances the 
XiXi^^ coupling, which is proportional to the difference of the squares of these components 

Figs. [2] show results analogous to those in Fig. [H except that now the gaugino mass 
parameter M1/2 has been varied along with Mx such that the LSP mass is kept fixed. This 
required taking larger values of M1/2 for smaller Mx- As expected from our previous discussion, 
the effect of reducing Mx is now quite small. The spin-independent cross section (left frame) 
increases by ~ 20% as Mx is reduced to its minimal value. This can be explained as follows. 
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Figure 2: As Fig. [H except that M1/2 has also been varied along with Mx, such that the mass 
of the lightest neutralino is kept fixed at 140 GeV. 



Since now m^o is kept fixed, we also have to keep fi essentially fixed in order to maintain the 
correct relic density. This requires reducing mo when Mx is reduced. This in turn leads to 
a reduction of rriA, which over-compensates the increase of that would result if Mx were 
reduced for fixed mo and fixed LSP mass. This implies a similar reduction for the mass of 
the heavier CP-even neutral Higgs boson whose exchange plays a prominent role in this cross 
section. However, it is not clear whether this variation is significant given astrophysical and 
likely experimental uncertainties. The variation of the spin-dependent cross section is even 
smaller. 

Figs. |3] shows the same cross sections for smaller (heaviest) neutrino mass, m^, = 0.2 eV, 
as well as larger gaugino mass, M1/2 = 1 TeV. Recall that the smaller m,^ requires a larger 
Yukawa coupling Yn, which, among other things, reduces the weak-scale f masses. This allows 
to satisfy the relic density constraint ([3]) for two distinct choices of mo- We continue to call 
the choice with ml ^ ^1/2' resulting sizable higgsino component of the LSP, the "focus 
point" [1^, even though the SO{10) model does not show "focusing" behavior of any Higgs 
soft breaking mass [3]. In the "co-annihilation" region the relic density is largely determined 
by Xi ~ ^1 co-annihilation [15] in both mSUGRA and the S'O(IO) model. 

For the focus point, we find the SI cross section to be almost independent of Mx- Note 
that m^o is now well above m^ in the entire range of Mx shown. Moreover, m^ increases with 
decreasing m^, the decrease of rriA with decreasing Mx is therefore less pronounced than 
in Fig. [H Finally, reducing the LSP mass increases the annihilation cross section, which scales 
like uiZq away from thresholds. In compensation, gaugino-higgsino mixing has to be reduced. 

Xi 

This reduces the LSP couplings to neutral Higgs bosons, offsetting the effect of the reduction 
of m^ as far as the SI cross section is concerned. In the SD case, we again observe a slight 
increase of the cross section with decreasing Mx, as in Fig. [H away from the tt threshold. 
Note also that, in spite of the increased LSP mass, the "focus point" scenario remains easily 
testable by near-future direct search experiments, at least via the SI cross section. 
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Figure 3: Dependence of the SI (left frame) and SD (right frame) neutrahno-proton scattering 
cross section on Mx, for rrii, = 0.2 eV. We chose two points that represent the "focus point" 
("FP") region (black) and the "co-annihilation" ("CO") region (red or grey). The neutralino 
mass varies between 220 GeV and 390 GeV at the smallest and largest Mx, respectively. The 
experimental constrains are taken for m^o = 220 GeV, and scale essentially like fn^o. 



On the other hand, for the co-annihilation point, both the SI and SD cross sections increase 
by one order of magnitude when Mx is reduced to its smallest allowed value. Here the Dark 
Matter relic density is mainly determined by the mass difference between the LSP and the 
lightest stau, which does not strongly depend on /i. Instead, the correct relic density is 
obtained through the direct effect of mo on m^^ . Due to the strong (exponential) dependence 
of the relic density on the Xi ~ mass splitting, only relatively minor adjustments of mg are 
required, which do not lead to significant changes of fi. In contrast, the additional Yukawa 
couplings in the 5'O(10) model reduce all over the parameter space. Therefore, in the 
co-annihilation region the larger higgsino component of Xi gives rise to larger scattering 
amplitudes, in particular via the Higgs- and exchange diagrams that dominate the SI and 
SD cross sections, respectively. The SI cross section is enhanced in addition by the decreasing 
rriA- As a result, at the smallest value of Mx this cross section even approaches the XENONIOO 
sensitivity. 

In order to understand the strong dependence of these cross sections on Mx, one has to 
keep in mind that reducing Mx increases the effect of the new Yukawa coupling in two ways. 
First, reducing Mx reduces the intermediate scales and Mq even more, i.e. \ia{Mx / Mc) 
and ln{Mx/Mji) increase when Mx is decreased. This increases the energy range where this 
coupling is effective in the RGE. Secondly, the reduction in Mji has to be compensated by 
the increase of in order to keep the very large Majorana neutrino mass in the see-saw 
expression ([T]) constant. 

The most robust indirect neutralino Dark Matter detection signal [10] is due to the capture 
of DM particles by the Sun, which greatly enhances the neutralino density near the center 
of the Sun. Eventually capture and annihilation of neutralinos in the Sun reach equilibrium. 
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Figure 4: Neutrino-induced muon flux from neutralino anniliilation in tlie Sun as a function 
of Mx for nil, = 0.4 eV (left); m,^ = 0.2 eV (right). The IceCUBE sensitivity limits are for the 
smallest LSP mass in the respective frames; in the relevant range of masses, the sensitivity 
limit scales roughly oc 1/m^o. 



The only annihilation products that can escape the Sun are neutrinos. In particular, muon 
neutrinos produce muons via charged current interactions; these muons can be searched for 
by "neutrino telescopes". 

In Figs. Hlwe plot the resulting muon flux as function of Mx, and compare it to the "best 
case" sensitivity of IceCUBE [16], using the input parameters of Fig. [1] (left frame) and Fig. |3] 
(right). Note that the overall neutrino flux is essentially fixed by the capture rate. The 
neutralinos interact with nuclei in the Sun mostly via Higgs and exchange. The capture 
rate is thus again sensitive to the higgsino components of the mostly bino-like neutralinos. It 
also depends on the mass of the neutralinos: the heavier the LSP, the less likely it is to lose 
enough energy in the interaction to become gravitationally bound to the Sun. The predicted 
muon flux therefore increases faster with decreasing Mx than the cross sections shown in 
Figs. [Hand [3] do. 

The muon flux also depends on the (mean) neutrino energy, since the neutrino charged 
current cross section increases with energy. Annihilation into pairs of or bosons leads 
to the hardest neutrino spectra, and hence to the largest signals. Annihilation into tt gives 
a somewhat softer spectrum, since some of the energy is taken away by the 6— quarks. This 
enhances the effect of the ti threshold visible in the left frame: To the left of this threshold, 
neutralinos predominantly annihilate directly into massive gauge bosons, while to the right of 
the threshold, annihilation into ti dominates. 

Of course, the neutrino energy also scales with the mass of the annihilating neutralinos. In- 
deed, the sensitivity limit on the muon flux decreases with increasing LSP mass for m^o < 500 
GeV [16]. However, in the muon flux itself this effect is compensated by the reduction of 
the neutralino flux impinging on the Sun, which scales like 1/m^o. Nevertheless, this effect 
keeps the expected flux in the "focus point" region well above the sensitivity limit even for 



the larger value of M1/2 chosen in the right frame. However, the flux in the co-annihilation 
region remains well below the IceCUBE sensitivity even for the smallest possible value of Mx- 
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Figure 5: SI neutralino-proton scattering cross section (left) and neutrino-induced muon flux 
from neutralino annihilation in the Sun (right), as function of m^. 



In Fig. |5l we show the SI proton-neutralino cross section as well as the neutrino-induced 
muon flux as function of mj,, for Mx = 3 ■ 10^^ GeV. We see that in the FP region, the 
cross section slightly increases with decreasing rrty. Recall that decreasing my, i.e. increasing 
the coupling Yat, reduces /x. In order to keep the relic density fixed one has to increase /i 
again by decreasing mo, which in turn leads to the decrease of the Higgs boson masses; this 
overcompensates the increase of tua with decreasing irii, if all soft breaking parameters are kept 
fixed. The reduced Higgs boson masses increase the scattering cross section. However, it also 
increases the importance of the A— exchange contribution to the annihilation cross section 
at rest. For tan/3 ^ 1, A— exchange mostly leads to hh final states, which produce very soft 
neutrinos. This effect over-compensates the (small) increase in the neutralino capture cross 
section, leading to the (very slight) decrease of the muon flux with decreasing mu in the FP 
region. 

In the co-annihilation region, increasing the Yukawa coupling Y)v reduces rrif^ as well 
as /i. The two effects tend to cancel, but a net reduction of nif^ results. This has to be 
compensated by increasing mo in order to keep the relic density in the desired range. This, as 
well as the effect of Yn in the RGE, increases itla- The increase of itia and the decrease of /i 
essentially cancel in the SI cross section. However, increasing itla also reduces the importance 
of neutralino annihilation to hb. This increases the average neutrino energy, which explains 
the slight increase of the muon flux with decreasing rriy. 

We have also computed the near-Earth flux of antiprotons due to the annihilation of relic 
neutralinos in the halo of our galaxy. As well known, the flux depends sensitively on several 
poorly known astrophysical quantities. One of these is the density of Dark Matter, which is 
reasonably well known "locally" , but not near the center of the galaxy, where it is largest. Note 
that, unlike positrons, antiprotons can diffuse from the galactic center to Earth. We illustrate 



this uncertainty by comparing three different halo models. The "N03" profile has been derived 
[T7] starting from a profile extrapolated from A^— body simulations [18], assuming that baryon 
infall compresses the Dark Matter distribution near the galactic center adiabatically. In the 
opposite extreme, one can assume that baryon infall heats the dark halo, leading [T7] to a 
profile similar to the (phenomenologically apparently quite successful) "Burkert" profile [T^ . 
Finally, the "NFW" profile [20] lies between these extremes. 

Figs. [6] show the dependence of the antiproton flux on Mx- Antiprotons are produced in 
the Galactic halo due to the hadronization of antiquarks produced in neutralino Dark Matter 
annihilation. As a result, the typical p energy is well below m^o. We show their differential flux 
at a kinetic energy 20 GeV, where the signal-to-background ratio is expected to be optimal 
[2T] . We illustrate the dependence on the halo model using the three profiles discussed above. 

The left frame of Fig. |6] is for the "focus point" region, with small Yn and relatively 
small Ml/2- In this case the relic density is determined by x? annihilation with itself, and is 
dominated by annihilation from the S*— wave, which is the only contribution relevant for the 
p flux. As a result, the x? annihilation cross section remains essentially constant in the left 
frame. However, we saw in Fig. [H where the same parameters were used, that m^o decreases 
by nearly a factor of two as Mx is decreased. This increases the x? annihilation rate, computed 
as the product of flux and cross section, by almost a factor of four. However, decreasing m^o 
also makes it increasingly more difficult to produce antiprotons at 20 GeV. As a result, the p 
flux near Earth only increases very slightly as Mx is decreased. 

The right frame shows results for a point in the co-annihilation region, with larger M1/2 
and smaller niy. Here the relic density is essentially determined by fi — x\ co-annihilation. 
The annihilation cross section increases significantly with decreasing Mx-, due to the decrease 
of (almost) all weak-scale sparticle and Higgs boson masses. Moreover, m^o now remains so 
high that getting 20 GeV antiprotons is not difficult. As a result, the rate increases by about 
an order of magnitude as Mx is reduced to its lower bound. 
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Figure 6: Antiproton differential flux for different halo models for = 0.4 eV ("focus-point", 
left) and rrii, = 0.2 eV ("co-annihilation", right). 




Figure 7: Regions in the (Mi/2,mo) plane for mu = 0.4 eV (left); rrii, = 0.2 eV (right). The 
grey area is excluded by the EWSB condition or by too light sfermions, and the scarlet area 
is excluded by the LEP limits on Higgs and chargino masses. The black points give us the 
correct Dark Matter relic density. The light blue region can be tested by searching for muon 
neutrinos originating from Xi annihilation in the Sun, while in the magenta region, the XiP 
scattering cross section exceeds the CDMS-II bound. The dependence of the detectability on 
the WIMP mass is taken into account, but we assume a fixed local WIMP density. 



This seems impressive, but is still smaller than the difference in the predictions based on 
the N03 and Burkert profiles. Additional systematic uncertainties come from the propagation 
of the antiprotons; here we have used DarkSUSY default parameters. Note finally that the 
p flux that can be inferred from the p/p ratio measured by the PAMELA satellite [22] and 
the well-known [5] proton flux is about 2 ■ 10~^ GeV~^cm~^s~^sr~^, well above even the most 
optimistic prediction in Fig. [61 Given that the prediction for the background also has sizable 
uncertainties, we conclude that the observation of cosmic antiprotons is not a very promising 
test of the models discussed here. 

Prospects for direct and indirect DM search in the S'O(IO) model with the smallest allowed 
Mx are summarized in Figs. [TJ for two different values of m^. The regions of parameter space 
that give the neutrino-induced muon flux from Xi annihilation in the Sun above the IceCUBE 
sensitivity limit are depicted as light blue. The regions where the spin-independent neutralino- 
proton cross section exceeds the CDMS-II bound are shown in magenta. Note that we always 
assume fixed local neutralino density when deriving these bounds, independent of the value of 
^yoh? predicted in standard cosmology. We also show the region excluded by the electroweak 
symmetry breaking (EWSB) condition or by too light sfermions (grey) as well as that excluded 
by the LEP limits [5j on the masses of Higgs bosons and charginos (scarlet). The black points 
are where the Dark Matter relic density satisfies Eq.(l3]). 

We find that, as in mSUGRA [23j, the region of high mo, where Xi has a sizable higgsino 
component, will soon be covered by direct searches and also by the neutrino indirect search. 
Recall from Figs. [5] that the region of parameter space to be probed by XENONIOO is much 



larger than that probed by IceCUBE. Compared to mSUGRA, for given M1/2, Aq and tan/3 
this region occurs at significantly lower values of mo; this is true in particular for small m,^, 
i.e. sizable Yukawa coupling Y/v (right frame). Moreover, in this region neutralino Dark 
Matter remains detectable out to much larger values of M1/2 than in mSUGRA, since the 
ratio m^o/Mi/2 is nearly two times smaller in our scenario than in mSUGRA. 

The co-annihilation region is difficult to see in Figs. [TJ since it is very narrow. It extends 
to Ml/2 — 750 (1400) GeV for m^, = 0.4 (0.2) eV. Unfortunately this region will not be tested 
by near-future Dark Matter search experiments. However, we saw in Fig. [3] that the XiP 
scattering cross section exceeds that in mSUGRA by about an order of magnitude. Much of 
this region will therefore be testable by ton-scale direct Dark Matter detection experiments. 

3 Collider searches 

In this section, we consider signatures at the LHC for our model. It remains sufficiently 
similar to mSUGRA that the overall search prospects, i.e. the reach in sparticle masses, is 
essentially the same in both models; i.e., discovery should be possible out to ~ 3 (2) TeV 
for rriq ~ {^)mg |23], once the LHC reaches its full energy and luminosity. Of course, at 
the smallest allowed value of Mx the reach in M1/2 is nearly two times larger in our model, 
but Ml/2 is not measurable by TeV-scale experiments. We will therefore focus on ways to 
distinguish the model from mSUGRA using measurements at colliders, with emphasis on the 
LHC. 

3.1 Benchmark points 

We performed detailed analyses of collider signals for two distinct benchmark points. The 
input parameters as well as superparticle and Higgs spectra are listed in Table [H We chose 
points that satisfy all constraints, including the Dark Matter relic density constraint (but 
ignoring the indication of a deviation of the magnetic dipole moment of the muon from the 
Standard Model prediction [25] )■ We chose Mx at its lower bound of 3 ■ 10^^ GeV, and small 
TUy = 0.2 eV, in order to maximize the differences between our model and mSUGRA. On the 
other hand, we chose the parameters of the mSUGRA points such that the sparticle spectra are 
as similar as possible to those of the corresponding 5*0(10) benchmark points. In particular, 
we adjust the values of M1/2 such that the gluino masses are essentially the same in both 
models. Moreover, we chose the same tuq in both models, since this gives similar first and 
second generation squark masses. In this way we hope to isolate the non-trivial effects of the 
additional couplings via the RGE. 

Point 1 is chosen such that, at least in the 5*0(10) model, the lightest neutralino has a 
significant higgsino component. This requires mo > M1/2 even in this model. However, for the 
same gluino mass, one would need much larger mo to achieve a similarly small fi in mSUGRA. 
This would put squarks out of the reach of the LHC, making the scenario easily distinguishable 
from our 5*0(10) point. We instead chose to increase tan/? from 40 to 52, and also took a 
nonvanishing (but fairly small) Aq. This leads to greatly reduced mass of the CP-odd Higgs 
boson, i.e. we are now close to the "A— pole" region [26j where Xi annihilation is enhanced 
since A— exchange in the s— channel becomes (nearly) resonant. These changes do not affect 
fi very much, i.e. in our mSUGRA point 1 the LSP remains a nearly pure bino. 

Benchmark point 2 lies in the co-annihilation region. Recall that the new coupling Y/v 
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Table 1: Benchmark points used in our analysis of collider signals. Mass spectra are calculated 
using a modified version of SOFTSUSY 2 . 0. Rows 2 through 5 give the input parameters. Row 
6 is the weak-scale value of Rows 7 through 28 give the on-shell mass of the indicated 
sparticle or Higgs boson, while row 29 gives the prediction for the scaled LSP relic density. 
Finally, the last row is the longitudinal polarization of the r lepton in fi t~x5 decays. All 
dimensionful quantities are in GeV. 

reduces m^^ below its mSUGRA prediction. Choosing M1/2 such that one gets the same g (or 
Xi) mass, while keeping all other input parameters the same, would thus lead to an mSUGRA 
point with too high a rehc density. We consider two different methods to correct for this. In 
mSUGRA point 2a we take non-vanishing Aq, such that m^^ is reduced and /i is increased; the 
latter also decreases m^^, helping to get a sufficiently large Xi~^i co-annihilation cross section. 
In mSUGRA point 2b, this is instead achieved by increasing tan/3, which again reduces mf^ 
and increases tl — tr mixing. Notice that in either case the change of these input parameters 



is not very dramatic. 

At the SO{10) point 1, all of the squarks as well as the gluino are significantly heavier than 
all of the neutralinos and charginos. Furthermore, due to the low value of |yu|, the heaviest 
neutralino Xa has the largest SU{2) gaugino component. Due to the small higgsino mass, 
the dominant SUSY production channel is qq — xtxf- Mostly due to this process, the total 
inclusive SUSY production cross section at ^/s = 14 TeV is nearly three times larger than 
that of the mSUGRA point 1. However, X23 xt decay predominantly into Xi ^ 
quark-ant iquark pair, which carries relatively little energy due to the small mass splitting. 
Direct X2 sxf production therefore predominantly gives rise to events with four relatively soft 
jets, and correspondingly only a small amount of missing Et- This signal will be completely 
swamped by backgrounds, e.g. from W, Z plus multi-jet production. The inclusive cross 
section for squark and gluino production, which should be detectable in this scenario (see 
below), is quite similar in the 5*0(10) and mSUGRA versions of point 1. 

Another distinctive feature of 5*0(10) point 1 is the much smaller polarization Pr of r 
leptons produced in ff — )■ r~Xi decays. Pr depends [27] both on tl — tr mixing and on 
gaugino-higgsino mixing. In the case at hand, fi is dominated by the tr component in both 
the 5*0(10) model and in mSUGRA; the tl component is slightly smaller in the 5*0(10) case 
due to the reduced value of /itan/3. However, the 5*0(10) model features much stronger bino- 
higgsino mixing in this case. Note that the bino couples to r^, while the (down-type) 
higgsino couples to tl. As a result, Pr is significantly smaller in the 5*0(10) case. 

Pr can be measured via the energies of hadronic r decay products [2S]. Of course, this 
requires a copious source of fi particles. At an e^e~ collider this measurement can therefore 
only be performed if the beam energy is well above m^^, i.e. ^/s > 2 TeV in our case. Monte 
Carlo simulations indicate [29] that Pr could then be determined with sufficient accuracy to 
distinguish these scenarios. At the LHC this measurement is probably only possible if fi 
particles are produced copiously in the decays of gluinos and/or squarks [30j; this is not the 
case in our benchmark point 1. 

For the 5*0(10) point 2, since the gluino and squarks are relatively lighter, the dominant 
sparticle production process is qg — )■ gqL,R- Our choices of M1/2 and mo ensure that the 
corresponding cross section is very similar in the 50(10) and both mSUGRA scenarios. 

Recall that we adjusted the mSUGRA parameters such that we get very similar m^^, and 
hence similar LSP relic density, as that in the 5*0(10) scenario. These adjustments also imply 
that the masses of third generation squarks are only slightly smaller in the 50(10) scenario 
than in both mSUGRA scenarios, i.e. the effect of the new Yukawa couplings on sfermion 
masses has been partly compensated by adjusting soft breaking parameters. However, the 
effect of the new couplings is still visible in which is significantly smaller in the 5*0(10) 
benchmark point than in both mSUGRA variants. 

Moreover, having adjusted parameters such that we obtain similar gaugino and first gen- 
eration squark masses, we get significantly heavier first generation sleptons in 5*0(10) than in 
mSUGRA This can be tested trivially at e'^e~ colliders operating at ^/s > Irrie^. However, 
even in the mSUGRA versions of our point 2, sleptons are too heavy for direct slepton pair 
production to yield a viable signal at the LHC [3T] . 

In the next Subsection, we will therefore focus on events containing charged lepton pairs 
originating from the decays of squarks and gluinos. We will show that this allows to distinguish 
the 5*0(10) points from their mSUGRA analogues, using the fact that is smaller in the 
5*0(10) scenarios. 



3.2 Measurements using di— lepton events at the LHC 



In order to analyze the gaugino-higgsino sector of the theory, we have to rely on neutralinos 
and charginos produced in the decays of squarks and gluinos. Direct production of charginos 
and neutralinos is only detectable in purely leptonic final states |32]. In the case at hand the 
relevant neutralino and chargino states are quite massive, and have small leptonic branching 
ratios, leading to very small signal rates. 

We therefore look for events with several energetic jets in addition to two or more leptons. 
To that end, we simulate proton-proton collisions at the LHC {y/s = 14 TeV) using PYTHIA 
6.4 [23! detector PYCELL. The detector is assumed to cover pseudorapidity \ri\ < 

5 with a uniform segmentation At] = A0 = 0.1. We ignore energy smearing, which should 
not be important for our analyses. We use a cone jet algorithm, requiring the total transverse 
energy Et summed over cells within R = 0.4 to exceed 10 GeV; here R = \/ [Sr])'^ + (50)^, 
where 61] and 6(f) measure the deviation in pseudorapidity and azimuthal angle from the jet 
axis, tt and diboson production are assumed to be the main Standard Model backgrounds in 
the di-lepton channels we are interested in. 

We require electrons and muons to have pt > 10 GeV and to be isolated, i.e. to have less 
than 10 GeV of additional Et in a cone with R = 0.2 around them. Also, leptons within 
R < 0.4 of a jet are not counted. These requirements essentially remove leptons from the 
decay of c and b quarks. Finally, we require the invariant mass of opposite-sign, like-flavor 
lepton pairs to exceed 20 GeV, in order to suppress contributions involving virtual photons. 

3.2.1 Point 1 

In this case the difference in fi between the 5*0(10) and mSUGRA scenarios is quite drastic: in 
the S'O(IO) case, is only slightly above Mi and well below M2, leading to m^o ~ m^o ~ m^±, 
only about 70 GeV above the LSP mass, and well below the masses of the wino-like xf and 
X4 states. In contrast, in the mSUGRA scenario we have fi slightly above M2, leading to 
wino-like X2 xf well above the LSP. 

In order to understand what this means for multi-jet plus di-lepton signatures, we have 
to analyze the most important sparticle production and decay channels. In the case at hand, 
the most important production channels (after cuts) are squark pair and associated squark 
plus gluino production, where the squarks are in the first generation. Most squarks will decay 
into a gluino and a quark here, so that most events start out as gluino pairs with one or two 
additional jets. 

In the 5*0(10) version of point 1, nearly all gluinos decay into ti plus top, since this is the 
only allowed two-body decay of the gluino. In turn, ii decays mostly into Xi^ and Xi 2 3^- 
These decays are preferred by phase space, and because here the /2(7/iter neutralinos and lighter 
chargino are dominantly higgsino-like, and hence couple more strongly to (s)top (since the top 
Yukawa coupling is larger than the electroweak gauge couplings). Leptons can then originate 
from semi-leptonic decays of top quarks, from leptonic decays of xf states, and from leptonic 
decays of X2 3- Note that the latter decays, which have branching ratios near 3%, can only 
produce di-lepton pairs with invariant mass below 70 GeV. 

In contrast, in the mSUGRA version of point 1, gluinos can only undergo three-body 
decays. Nevertheless decays involving third generation quarks are strongly preferred, since 
the bi and ti exchanged in g decay can be nearly on-shell. The dominant decay modes 
again involve higgsino-like states, i.e. g xt^i or x^itt] due to the larger phase space, the 



branching ratio for g — > Xibi is also significant. Tlie liiggsino-like states decay into lighter 
gaugino-like states plus a real gauge or Higgs boson. Leptons can then originate from semi- 
leptonic top decays, and from the decays of the and Z° decays produced in the decays of 
the heavier neutralinos and both charginos. Note that we do not expect any structure in the 
di-lepton invariant mass below Mz in this case. 

We apply the following cuts to suppress the Standard Model background [M] : 

• At least four jets with Et > 150 GeV each, at least one of which satisfies Et > 300 GeV. 

• Missing Et > 200 GeV. 

• Transverse sphericity St > 0.2. 

• Two charged leptons with opposite sign and same flavor (OSSF). 

No SM diboson event passed these cuts, and one ti event passed, for a simulated integrated 
luminosity of 1 fb~^. Note that our cuts are quite generic, not optimized for our scenario. In 
our case, the background can be further suppressed by requiring at least three tagged b quarks 
in the event (all the final states we discussed above have at least four b quarks); by requiring 
the presence of additional jets (most events will have at least one hard jet in addition to the 
gluino pair, which by itself already produces at least four jets); and/or by optimizing the 
numerical values of the cuts employed. This should allow to extract an almost pure SUSY 
sample, without significant loss of signal. 
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Figure 8: Left frame: The dilepton invariant mass distribution of opposite sign, same flavor 
dilepton events after cuts for point 1, corresponding to an integrated luminosity of 300 fb~^; 
the black and red (grey) histograms are for the mSUGRA and S'O(IO) versions of this point, 
respectively. Right Frame: The subtracted (OSSF— OSOF) dilepton invariant mass distribu- 
tion after cuts for point 2; the solid black and dashed blue histograms are for the mSUGRA 
points 2a and 2b, while the solid red (grey) histogram is for the SO{10) version. 



The results of our simulations for point 1 are shown in the left frame of Figs. |H1 We see 
that the di-lepton invariant mass distribution peaks near 50 GeV in the 5*0(10) scenario. 



whereas it peaks at Mz in the mSUGRA scenario. mSUGRA also predicts a somewhat larger 
number of events at large di-lepton invariant mass; we saw above that only in this scenario 
we expect significant numbers of on-shell bosons from chargino and neutralino decay. 
The two distributions should be easily distinguishable, with high statistical significance, once 
several hundred fb~^ of data will have been collected. 

3.2.2 Point 2 

We now turn to benchmark point 2. We saw in Table [1] that now X2 ~^ ^^^^ allowed in the 
mSUGRA scenarios, but not in S'O(IO). Unfortunately, the branching ratios for these decays 
remain at the permille level even in the mSUGRA scenarios. This is partly due to the small 
phase space available for these decays, but mostly due to the fact that X2 is dominantly a 
neutral wino in the co-annihilation region, and thus has only small couplings to SU{2) singlet 
sleptons; recall that fi also has a significant SU{2) doublet, tl component. As a result, we 
do not see any evidence for 

X2 ^^^^ ^ X?^^^^ = e, /i) in the mSUGRA scenarios; m 
particular, no kinematic edge at m£+£- = 189 GeV, the nominal endpoint for this decay chain, 
is visible. 

We instead first try to find the evidence for increased gaugino-higgsino mixing in the 
SO{10) scenario, due to the lower value of n, by analyzing the decays of the heavier charginos 
and neutralinos. In the case at hand the dominant production channel is associate production 
of a first generation squark with a gluino. The decays of first and second generation squarks 
will predominantly produce Xi 2 xf states. On the other hand, here in the co-annihilation 
region gluinos are heavier than all squarks, but - as usual in scenarios where squark masses 
unify at some high scale [33] - gluino decays into third generation quarks and squarks are 
preferred. We saw in the discussion of benchmark point 1 that third generation squarks in 
turn frequently decay into higgsino-like charginos and neutralinos. We look for these heavier 
states through their decays into real bosons. 

The relevant branching ratios are summarized in Table |2i We see that the slightly reduced 
third generation squark masses of the 5*0(10) scenario increase the branching ratio for gluino 
decays into the third generation to 69% in the 5*0(10) case, compared to 61% (60%) in 
mSUGRA point 2a (2b). Moreover, the decays of third generation squarks into Xs ^ind xt 
are enhanced in the 50(10) case. This is also predominantly a phase space effect; the slightly 
reduced squark masses in the 50(10) case are over-compensated by the reduced masses of 
the higgsino-like states. This effect is especially drastic for ti — )■ Xgt, where the available 
phase space volume is very small in the mSUGRA scenarios. However, the strong phase space 
dependence of the relevant partial width^ [6j leads to quite significant differences also in 
many other modesfl Finally, the smaller value of /i in the 50(10) scenario also leads to more 
higgsino-gaugino mixing, and hence to slightly larger branching ratios for decays of the heavier 
neutralinos and charginos into real bosons. The combination of these three effects leads 
to a substantially larger Z° production rate in gluino cascade decays in the 50(10) scenario 
than in mSUGRA. 

*If niq is negligible, the partial width for g — ^ x + 9 is oc (1 — m?/™?)^. 

^This also holds for decays into X4- However, Br(x4 — >■ Xi 2-^°) '^^^Y amounts to ^ 3%. This large difference 
between the decay modes of the higgsino-like states X3 and Xa can be traced back to the fact that x-j is a very 
pure symmetric higgsino, i.e. the higgsino components of this eigenvector are nearly equal in both magnitude 
and sign, whereas xi is mostly an antisymmetric higgsino. As a result, the XaXi 2-^*' couplings are suppressed, 
and the X4X1 2^ couplings are enhanced, where h is the light neutral Higgs boson. 
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6.2 % 
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8.9 % 
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Table 2: Branching ratios for important modes in benchmark point 2, as calculated with 
ISA JET 7.78 [36j; note that charge conjugate gluino decay modes have to be added. The last 
line denotes the sum of the branching ratios of all gluino decay chains which give us a Z boson 
in the final state. 

In order to suppress backgrounds, we first make use of the fact that most signal events 
will have (at least) one very energetic jet from the decay of a first generation squark into a 
light gaugino-like state. Moreover, the above discussion shows that many events with real Z^ 
bosons in the final state also will have a top quark in the final state, and/or a real boson 
from xt decay; the branching ratio for this latter decay amounts to 12.6% in the S'O(IO) 
scenario. The decays of these particles can lead to additional, somewhat softer, jets and/or 
additional leptons. On the other hand, since the dominant production channel only contains 
a single gluino, we only expect two b (anti-)quarks in the final state; 6— tagging will therefore 
not be of much help to suppress the tt background, which we again expect to be the most 
dangerous one. 

This leads us to use two complementary sets of cuts; at the end we simply add both event 
samples in order to increase the statistics: 

1) Set 1 

• EtUi) > 600 GeV, ^t(J2) > 200 GeV, EtUs) > 100 GeV. 

• Ni = 2. 

2) Set 2 

• EtUi) > 300 GeV, Et(J2) > 150 GeV, EtUs) > 75 GeV. 

• > 3. 

Note that we do not apply any cut on missing Et, since this wasn't necessary to suppress 
the backgrounds we studied. In Set 1, a modest missing Et cut will be needed to suppress 
the Z°+jets background, but this can be done without any significant loss of signal. 



For set 1, we require both charged leptons to be of opposite charge. In set 2, we chose the 
opposite-sign lepton pair whose invariant mass is closer to Mz as candidate" . We found 
that the background from tt is almost entirely removed by either set of cuts. Furthermore, in 
order to extract the lepton pairs from the decays of Z° bosons, we subtract the events with 
opposite sign opposite flavor (OSOF) lepton pairs. This removes SUSY backgrounds where 
the two leptons originate from independent (semi-)leptonic decays; in benchmark point 2, 
these come primarily from the decays of bosons. 

The resulting di-lepton invariant mass spectrum is shown in the right frame of Fig. [HI for an 
integrated luminosity of 300 fb~^. We see that the peak is indeed much more pronounced in 
the SO{10) scenario than in the two mSUGRA scenarios. This allows to distinguish between 
5*0(10) and mSUGRA at about 3cr statistical significance in this case. 

As noted above, in the 50(10) point 2, we find a branching ratio for — )■ W^Xi of 
about 12.6%. In the mSUGRA points 2a and 2b, this branching ratio is only 6.1% and 6.6%, 
respectively. At an e+e~ linear collider with sufficient energy to produce XiXi pairs this large 
difference in branching ratios should be straightforward to measure. 

At the LHC we have to pursue a somewhat different strategy: leptonic decays of these 
can give rise to events with two leptons of the same charge (like-sign di-lepton events). 
These can originate from associate qig production where g decay also produces a lepton; the 
charge of this lepton from gluino decay is uncorrelated to that from squark decay, i.e. half 
the time the two leptons will have the same charge. The results of Table [2] indicate that the 
inclusive branching ratio for g ^ is somewhat higher in the 50(10) scenario than in both 
mSUGRA analogues. Other sources are ulUl and didi production, which give rise to 
and pairs, respectively, if both squarks decay into x^- Of course, gluino pairs can also 
produce like-sign dileptons, but the gluino pair production cross section is relatively small at 
this benchmark point. Note also that the physics background for like-sign dileptons events is 
very small. 

We applied the following cuts to isolate a clean sample of SUSY events: 

• At least two jets, with Et^Ji) > 500 GeV, -E'T(j2) > 200 GeV. These cuts are quite 
asymmetric, since we expect at least one very energetic jet from decay in the event. 

• Exactly two equally charged isolated leptons, with pt{£) > 10 GeV as before. 

With these cuts, we find 492 events in 300 fb~^ for the 50(10) scenario, as compared to 422 
and 434 events in mSUGRA 2a and 2b, respectively. This difference of ~ 3 statistical standard 
deviations is much less than the above discussion would lead one to expect. The reason is that 
there is another large source of £^ from decay: — )■ gx2 with X2 ~^ t'^'^i^, and — )■ (.^uu^ 
Unfortunately the branching ratios for these decays are somewhat larger in mSUGRA than 
in the 50(10) scenario, because higgsino-gaugino mixing tends to suppress the corresponding 
partial widths. This significantly reduces the difference between the predictions for the total 
like-sign dilepton event rate. 

One can imagine two strategies to enhance the difference between the predictions. One 
possibility is to veto leptons from X2 ~^ t^t^ decays by vetoing against the secondary 
from fi decay. However, this r will be quite soft, so it is not clear how efficient such a r 
veto would be. Another possibility is to subtract this source of hard leptons, by using events 
with an identified r jet and the known r decay branching ratios. Again, the feasibility of this 

■'"The decay of the fi produced in this chain, or via the dominant decay — >■ f^v, only produces very 
soft T leptons, and hence even softer since we are in the co-annihilation region where the — x\ mass 
difference is small. 



method depends on r tagging efficiencies and their uncertainties. We do therefore not pursue 
this strategy any further. 

3.3 Measurements involving Higgs bosons at the LHC 

Our benchmark points have quite large values of tan /3. This increases the cross sections for 
inclusive gg A, H production, and for associate gg — bb{A, H) production. The heavy 
Higgs bosons can be searched for using their decays into r+r^. According to simulations by the 
CMS collaboration [37], for tan/3 = 40 this would allow discovery of the heavy Higgs bosons 
out to niA — 650 GeV with 60 fb^^ of data. In particular, we expect a robust signal for H, A 
production in the mSUGRA version of benchmark point 1, but not in the 5*0(10) version. In 
benchmark point 2, we expect signals of comparable size in all three cases. The t~^t~ invariant 
mass resolution should suffice to distinguish between the 50(10) and mSUGRA 2a scenarios, 
but distinguishing the 50(10) scenario from mSUGRA 2b might be challenging. 

In benchmark point 2, X2 ~^ + ^ decays might also allow to discriminate between 50(10) 
and the two mSUGRA analogues. The branching ratio for this decay is 11.5% in the 50(10) 
case, but only 5.3% (5.7%) for mSUGRA 2a (2b). About 90% of the light Higgs bosons will 
decay into bb pairs. Recall, however, that in this scenario gluino decays frequently lead to 
bb pairs in the final state, giving rise to a large SUSY background. We have therefore not 
pursued this avenue further. 

4 Summary and Conclusions 

We have investigated the detectability of neutralino Dark Matter and some LHC signatures for 
a SUSY— 5*0(10) model with two intermediate scales, and compared results with the frequently 
analyzed mSUGRA model. 

In Sec. 2 we have shown that in the cosmologically allowed region with large scalar mass 
parameter mo, the direct detection of Dark Matter should be possible for the next generation 
of detectors, at least for gluino masses up to 2 TeV. In this region of parameter space, which 
corresponds to the "focus point" region of mSUGRA, IceCUBE will be able to detect the 
neutrino-induced muon flux from neutralino annihilation in the Sun. However, in this region 
of parameter space the Dark Matter detection rates predicted in the 5*0(10) scenario are quite 
similar to those in mSUGRA, if one adjusts parameters such that the physical LSP mass and 
relic density is the same in both cases. On the other hand, in the co-annihilation region 
prospects for direct Dark Matter detection are much better in the 5*0(10) case, although 
the cross section still remains somewhat below the projected sensitivity of next generation 
experiments. Unfortunately, in this case the neutrino signal from the Sun is several orders of 
magnitude below the IceCUBE sensitivity. 

We also analyzed the flux of antiprotons from neutralino annihilation in the halo of our 
galaxy. Even in the most optimistic case - with a large annihilation cross section, and a 
halo proflle strongly peaked at the galactic center - the predicted flux falls well below the 
flux extracted from the PAMELA p/p flux ratio and the known proton flux. Given the large 
uncertainties in both signal and background, searches for antimatter do not appear to be very 
promising Dark Matter search channels in this scenario. 

In Sec. 3 we compared some collider signals from the model with those of mSUGRA. We 
chose two benchmark points, one of which resembles a "focus point" scenario, while the other 



lies in the co-annihilation region. We compared those with mSUGRA points with identical 
gluino and first generation squark masses, and (almost) identical neutralino relic density. We 
showed that events containing hard jets and two (or more) charged leptons (electrons and 
muons) can be used to discriminate between mSUGRA and our 5*0(10) benchmark points. 
The existence and size of a Z° peak in the di-lepton invariant mass distribution proved par- 
ticularly useful. Moreover, in the first benchmark point the search for heavy neutral MSSM 
Higgs bosons can help to discriminate between the SO{10) and mSUGRA models, whereas in 
the second point the number of like-sign dilepton events can be used. 

It may be interesting to note that searching for events containing leptonically decaying Z° 
bosons [2S] or like-sign dileptons ^39j were the first strategies suggested to look for cascade 
decays of heavy squarks and gluinos. It was realized later that other channels with fewer 
leptons offer better SUSY discovery reach [lOl [2l] . Here we find that bosons and like-sign 
lepton pairs can be very useful for discriminating between different SUSY models. 

In the co-annihilation region, the physical spectra predicted by 5*0(10) and mSUGRA are 
quite similar. The models differ mostly in the values of fi derived via electroweak symmetry 
breaking; the 50(10) scenario also predicts somewhat reduced masses for third generation 
squarks. The fact that our analyses nevertheless led to significant differences in several signals 
bodes well for the power of LHC experiments to distinguish between competing SUSY models. 
However, we have not attempted to distinguish the 50(10) model from other modifications 
of mSUGRA, e.g. scenarios with non-universal soft masses for Higgs bosons [41j, which also 
lead to variations in /i as well as the masses of third generation sfermions and Higgs bosons. 

Indeed, it is clearly impossible to distinguish our model from sufficiently complicated 
SUGRA scenarios with non-universal scalar soft breaking masses using superparticle and 
Higgs boson properties only, simply because the two models can have identical weak-scale 
spectra. Recall, however, that our model also makes other predictions. For the chosen (ex- 
treme) value of the scale of Grand Unification, proton decay should be within reach of next 
generation experiments [12]. This low unification scale also implies a rather large mass for at 
least one neutrino; this can be probed using cosmological data [33] , and perhaps even through 
laboratory experiments [H]. By combining all the information, we will hopefully be able to 
pin down the physics at the Grand Unified scale. 
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